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Abstract. The joint effect of the electron-phonon interaction and Kondo effect on 
the nonequilibrium transport through the single molecule transistor is investigated by 
using the improved canonical transformation scheme and extended equation of motion 
approach. Two types of Kondo phonon-satellites with different asymmetric shapes are 
fully confirmed in the spectral function, and are related to the electron spin singlet or 
hole spin singlet, respectively. Moreover, when a moderate Zceman splitting is caused 
by a local magnetic field, the Kondo satellites in the spin resolved spectral function are 
found disappeared on one side of the main peak, which is opposite for different spin 
component. All these peculiar signatures that manifest themselves in the nonlinear 
differential conductance, are explained with a clear physics picture. 
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1. Introduction 

The Kondo effect as a manifestation of strongly correlations between conduction 
electrons and spin impurity has been extensively studied in the context of quantum 
dot physics in past years. [1-7] Recently, its realization in the single-molecule transistor 
(SMT) has attracted a lot of attentions. [8-14] In particular, the latest experiments in 
SMT, in which the phonon satellite has been clearly observed even in the Kondo regime 
[15-17] owing to the enhanced electron-vibration coupling in the molecule, [18-20] has 
stimulated great interests to investigate the interplay of the electron-phonon interaction 
(EPI) and the Kondo effect on the transport through the SMT. 

According to our recent investigation on the spectral function of the SMT in the 
absence of the Coulomb interaction, the phonon satellites are quite sensitive to the lead 
Fermi levels with respect to the localized level. [21] When the temperature is too low 
to excite phonons thermally, the phonon satellites can only develop below (above) the 
resonant level if this level is occupied (empty) initially, as then electrons can only emit 
phonons. On the other hand, in the Kondo regime with the EPI ignored, it is well 
known that two collective spin singlet states will be formed between the localized state 
and the excited electron states above the Fermi level fx or the excited hole states below 
fx because of spin exchange processes, [22,23] which will contribute to the formation of 
the sharp Kondo peak pinning at the Fermi surface in the local density of states. It is 
thus expected that the interesting transport physics may result from the interplay of 
both the EPI and the Kondo effect, in particular out of the equilibrium, as both effects 
are sensitive to the Fermi levels at the source and drain electrode. 

Several theoretical techniques have been developed and applied to address the 
joint effects of the EPI and the Kondo correlation on transport through a quantum 
dot or a SMT. [24-34] According to these works, the equilibrium properties, in 
particular the renormalized effect on the Kondo correlation due to the EPI, have been 
well described. For example, with the help of the numerical renormalization group 
approach, the renormalized effect has been predicted for wide range of parameters in 
the equilibrium. [26-29] With the help of generalized Schrieffer- Wolff transformation, 
under the assumption that the system stays equilibrium in the strongly asymmetric 
coupling case, the nonlinear differential conductance has been given in Ref. [30] with 
the Kondo satellite structure exhibited. On the other hand, owing to the lack of 
satisfactory treatment for the nonequilibrium Kondo problem, little effort has been 
made to investigate how the Kondo phonon-satellites develop and manifest themselves 
in the nonequilibrium transport. However, since the Kondo phonon-satellites can only 
be observed at finite bias voltage experimentally, when the system is driven out of 
equilibrium, [15-17] it is certainly of importance to investigate the key signature of the 
nonequilibrium phono n-Kondo effect. As far as we know, by the real-time diagrammatic 
formulation, Koning et al. have given the spectral function as well as the differential 
conductance with Kondo satellites in the nonequilibrium situation. [24, 25] Here we 
shall develop another easier and straightforward approach, which will combine the 
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improved nonperturbative canonical transformation treatment of the EPI [21] with the 
extended equation-of-motion (EOM) method of the nonequilibrium Green functions. 
This approach, although is not rigorous in quantitative description, has the advantage 
of intuitiveness and can provide a semi- quantitative understanding of the phonon- 
assisted Kondo effect, particularly for strong electron- vibration coupling as well as 
nonequilibrium situation. 

Based on the approach mentioned above, in the present paper, we shall mainly focus 
on how the Kondo satellites develop and manifest themselves in the nonequilibrium 
transport through the SMT, and how these Kondo satellites are affected by an applied 
local magnetic field. After the Introduction we shall give our model and method in 
Section II and III, respectively. Our main results are presented in Section IV, which 
include: (i) The Kondo phonon-satellites may exhibit quite different asymmetric line- 
shapes with respect to the Kondo main peak in the spectral function, and also in the 
nonlinear differential conductance; (ii) Two types of spin exchange processes, associated 
with the collective spin singlets formed by electron states and by hole states respectively, 
can be clearly distinguished by the Kondo satellites; and (iii) With a moderate Zeeman 
splitting exceeded the width of Kondo peak, its Kondo-phonon satellites only appear 
on one side of the main Kondo peak in one spin-resolved spectral function, and on the 
opposite side for the opposite spin component. Finally, the conclusions are drawn. 

2. Model 

The SMT system studied in the present paper can usually be described as the Anderson- 
Holstein model, in which a single localized state coupled linearly to one local vibration 
mode and to the left (L) and right (R) non-interacting leads, namely 

H = ^2 £ akC^ kcT c ak(T + ^2 £ ° n ° + U^^i + ^oa f a 

a;k,cT u 

+ A M &t + a ) + (^<* c La + h - c -), (!) 

u a;k,cr 

where c^, ko .(c akcr ) and dJ.(d CT ) are the creation (annihilation) operators for the lead 
electron with energy e ak and the localized SMT electron with energy e a) respectively, a 
denotes the spin index and the lead index a = L,R. Uq is the on-site Coulomb repulsion, 
and n CT = dJ.d CT . The operator a^ (a) creates (annihilates) the local vibration mode with 
frequency c^o, A is the EPI strength, V^, k is the tunneling coupling between the localized 
and lead electrons, which results in a level broadening V = (T^ + r^)/2, where 

r a H = 27r^|K k | 2 ^-5 ak ). (2) 

k 

To treat the EPI non-perturbatively, by the canonical transformation, S = 
~ij-J2cr n °- ( at ~~ a )' the Hamiltonian is transformed into H = e s He _s = H p /i + H ei , 
where the phonon part tl p h = fkoos^a, and the electron part turns out to be the 
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Anderson Hamiltonian, namely 

a,k,cr a 

+ E(^c t Qk(7 d <T + /i.c.). (3) 

The parameters with bar in Eq.@ correspond to the renormalized ones resulted from 
the canonical transformation, e.g., the SMT level is shifted to e a = e a — gu>o with 
g = (X/ujq) 2 , and the on-site repulsion is renormalized to Uq = U — 2guj . When 
the charging energy is significantly reduced by screening due to the electrodes, [35] Uo 
may become negative for very strong EPI, which will result in the anisotropic Kondo 
effect. [26-29] But in most realistic situations, the charging energy is much larger than 
the EPI energy, [18-20] so we shall limit ourselves to the standard Kondo effect, i.e. TJq 
remains positive. Besides, = V^X. with X = exp [— (X/uq) (a* — a)], indicating 
that a cloud of phonons will be created or destroyed during the hopping processes. As 
containing the phonon operator X, the renormalized Hamiltonian H e ; has not been 
fully decoupled. For the present work, we are specially interested in the cases where the 
EPI is sufficiently strong compared to the tunneling coupling so that a local polaron is 
formed. Then we can approximate the phonon operator X by its expectation value in the 
thermal equilibrium, (X) = exp \—g (N p h + 1/2)], with N p h defined as the population of 
phonons at temperature T. Although it belongs to the mean field approximation, this 
method has been used widely and obtained reasonable successes. [36-41] 

With this mean field approximations, the lesser Green function can then be 
separated into the electron part and phonon part, 

G< (t) = i (4 (0) d CT (t)) = i (dle m d a e- m ) 
= i ^4e iHei *d CT e- iHei ^ ^ ^XV^Xe - ^**^ ^ . (4) 

The trace of the phonon part can be evaluated by the Feynman disentagling technique. 
[39] The lesser Green function of the SMT electron can be decoupled as 

oo 

G<H= £ UG^uj + nuo), (5) 

n=— oo 

and similarly the greater Green function is formulated as 

oo 

G>(u)= ]T L n G>(uj-nuj ). (6) 

71= — OO 

Here Ga is the dressed Green function associated with H e ;, and 

Ln = e -^N ph+ l) e n. Q ,2k B T In Ug^j Nph ( Nph + X )\ } ( 7 ) 

in which I n (z) is the n-th Bessel function of complex argument. The spectral function 
is then calculated with the lesser and greater Green functions via 

A a (u)=i(G>(u;)-G<(u J )). (8) 
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It should be emphasized that people usually evaluate the spectral function via A a {uj) = 
i{G r a {uj) —G"(u>)), in which the retarded and advanced Green functions, G r a and G°, are 
approximately separated into the electron and phonon parts by ignoring the difference 
between the N p h and N p h + 1, which works only in the high-temperature limit. [36-38] 
However, the phonon-Kondo problem studied here is certainly not the CclSG, clS the 
characteristic Kondo temperature, T K , is very low. Hence the prescription we have 
recently made may be more appropriate [21]. 



3. Keldysh Green functions 

Within the Keldysh formalism, the dressed Green functions G&^ can be derived from 
H e £ with the help of the EOM approach. Since the usual approximation in decoupling 
the higher-order Green functions generated by EOM is valid only nearby or above the 
Kondo temperature Tk, [42,43] at lower temperature as in the present situation the 
approximation developed by Lacroix may be more appropriate [44]. In the original 
Lacroix approach, only the system with the spin degeneracy and infinite charging energy 
was discussed in the equilibrium. To apply to the nonequilibrium situation with the finite 
Uq and Zeeman splitting, we have extended it accordingly. Thus, the usually ignored 
higher-order retarded Green function are taken into reconsideration. For example, 
((d^c^-c^k^d*)) is now approximately decoupled as (d 5 c^ k -) ((c a / kV , d+)), and the 
matrix element (d^c akd ) should be determined self-consistently by the nonequilibrium 
Green function technique. 

From the standard derivation of the EOM, [45, 46] the dressed retarded Green 
function turns out to be 

u - e a - £( T )(cu) + - (3) 

+ p-M) + p-A-" u ) (Q , 

u-e a -U - ggf H _ (3) 



where 



loI =uj + e<j- e a , 



u u = u - e a - e 9 - U Q , (10) 
and the correlation function between the SMT level and leads, P 9 (u>), is defined as 

p 9 (u) = ^t/* k (dk^>^ k H, (ii) 

ak 

in which g^i 00 ) * s t ne unperturbed retarded Green function of the lead electron. In the 
above the self-energy associated with the resonant tunneling processes is defined as 

& T \u;) = J2\^\ 2 9 r aM, (12) 

ok 
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and the self-energies related to the Kondo correlation are respectively expressed as 

= - s (T) M [p-M) + p;{-^)] , 



Here Y>f (uu) 



and 



F-Ju) 



£l 3) M 



0^ 



E tCkK'k'(cU'a C ak*R 



ok 



(13) 
(14) 



ak,o'k' 




Ba(w) = -2ilm { V* ly ., ( d^c a / k / s 



ak 



(15) 



In the wide-band limit, £( T )i 



£i 0) M 



2 i2vrA; B r 



E 



iY, and 

1 U U + fig 

2 i27rA; B T 



— 27T 



(16) 



where the Psi function ip is the logarithmic derivative of the gamma function. Now the 
self-consistent equation for Paiyj) is simplified as 

duo' [Y L f L {u')+Y R f R {u')]G^{u') 



P-M 



(17) 



2tt uo-oo' + i0+ 

which, together with Eq.(J£JJ), gives the dressed retarded Green function G r a . 

As we will see that at zero temperature, not only F s {u), but also Pa(uo) will diverge 
logarithmically and then all contribute to the sharp Kondo resonance at each Fermi 
surface. Therefore, Pa(co) should not be ignored as in the previous approaches. [42,43] 
Generally, this EOM approach works for most cases, however it fails to give the Kondo 
resonance at the particle- hole symmetric case. For example, when e-j- = fj_ = — Uq/2, 



CO 



v 



PZ(-uS) 



-P 9 (u), and F*(-uo) 



Therefore the divergent terms in the expression of the retarded Green function are 
cancelled with each other and then no Kondo resonance will show up. [47] The reason 
for this limitation is the fact that at the particle-hole symmetry point, the charge 
fluctuations are completely quenched, therefore the EOM approach can not be applied. 
[14] 

By the Keldysh formula Ga^ = G^Ha G%, in which the the greater (lesser) 
self-energy is evaluated by the ansatz adopted by Ng [48] , the dressed greater and lesser 
Green functions are obtained as 



G<(co) 



i([Y L (l - f L (co)) + Y R (1 - f R (co))})ImG r a (co)/Y, 
- i[Y L f L {u) + Y R f R {u)]ImG r a {u)/Y. 



(18) 
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Substituting Eq.(JTHJ) into Eqs.© and JBJ), one obtains the greater (lesser) Green function 
and subsequently the spectral function of the SMT. 

It is noted that, with the infinity charging energy and zero Zeeman splitting in 
the equilibrium, our results will recover to the Lacroix's one. [44] Together with the 
improved canonical transformation, this generalized EOM method will be our framework 
for investigating the Kondo satellites at the nonequilibrium with finite Zeeman splitting. 
Although this approach can only give a qualitative description of the shape of the Kondo 
satellites, it does predict the right peak positions. [26,42,43] For quantitative description 
of the renormalized properties in the equilibrium or exploring wider parameter space, 
more accurate methods like the NRG or quantum monte carlo calculations may be 
appropriate [26-29]. 

With the spectral function and the lesser Green function given above, the current 
through the SMT can be calculated via [49,50] 

J = E ^ / duJ M^ TL - fR^nlM^) 

+(T L -T R )G<(co)}, (19) 
and the differential conductance can be obtained straightforward as G = dJjdVuas- 

4. Spectral Function and Differential Conductance 

4-1. spin degenerate case 

First, we will discuss the spin degenerate case, i.e., £| = £j_ = s . The spectral function 
in the absence of the EPI is plotted in Fig^a) for different Uq, which shows that besides 
the two broad resonant peaks at Eq and Eq + Uq, there will be a sharp Kondo resonance 
developed at Hl(r)- By increasing Uq, as can be observed from the inset of Fig|T]^a), this 
Kondo resonant peak can be tuned more and more asymmetric, which is a manifestation 
of the particle-hole asymmetry of the SMT. 

When the electron-phonon interaction is turned on, as shown in Fig^b), the Kondo 
peak will be split into a set of phonon satellites spaced by multiple of phonon frequency. 
When the system is driven out of equilibrium by applied a bias voltage, the main 
Kondo peak will be split into two weakened peaks pinning at Fermi levels of the left 
and right lead, respectively, and the phonon Kondo satellites will also be split into two 
corresponding sub-sets. These results well agree with previous works, which validates 
the present framework [24,25,42,43]. 

The underlying physics about how and why the Kondo satellites develop in such a 
way can be intuitively understood with the help of Fig|5| As illustrated by FigE^al), a 
local hole state in the SMT exchanges with the excited holes at leads below the Fermi 
energy \x to form a Kondo spin-singlet. When an electron hops to the local state by 
emitting a phonon, as shown in FigEfa2), only the holes with energy below // — fkjQ 
can exchange with the local hole, because the phonon absorption is unavailable at very 
low temperature and the energy conservation law must be satisfied in the interaction 
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Figure 1. (color online) The spectral functions of the spin degenerate local electrons, 
i.e., A a , with (a) and without (b) the EPI. The parameters used in (a) are: hl — 
i = 0.0, £o = ~ 2I\ and Uq = 8T, Vbi as = 0. The parameters used in (b) are 
(/•»£ + Mfl)/2 = 0, £ = -2.5T, k B T = 2 x 10- 4 r, U = 100r, hu> = OAT. Other 
parameters have been shown explicitly in each figures. 




Figure 2. (color online) Schematic illustration for the spin exchange processes 
associated with the formation of the Kondo satellites. 
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Figure 3. The differential conductance for spin degenerate case as a function of the 
bias voltage. Here e = -2.5T, g = 0.5, k B T = 2 x lCT 4 r, U = WOT, and Huj = OAT. 
As an illustration, the inset shows the origin of the first Kondo phonon-satellite peak 
schematically. 



process. Thus a phonon-satellite peak at /x — Hujq appears in the density of states or 
tunneling spectra. Similarly, a spin-singlet consisting of a localized electron and excited 
electrons above /x is depicted in FigJ2fbl), and only the electrons above /x + fkv can 
exchange with the local electron if the hopping processes are accompanied by emitting 
a phonon as shown in Fig|2fb2), which will give rise to the satellite peak at /x + ftwo- 
That means although both processes shown in FigEfal) and Fig|2fbl) contribute to 
the main Kondo peak, they can be distinguished by their Kondo satellites, namely, the 
contributions of each type of exchange processes to the Kondo peak can be picked out 
in principle by comparing with the satellites on each sides of the main peak. 

Fig^b) also shows the satellites on different sides of the Kondo main peak may 
take on apparent sharp asymmetric lineshapes. The reason for these sharp feature lies 
in the sharp Fermi distributions in the leads at very low temperature, where the electron 
excitations and the hole excitations are separated from each other sharply by the Fermi 
surface. For example, the phonon-assisted spin exchange processes mentioned above 
have the sharp energy boundaries as Hl{r) ±nu , therefore they may result in the sharp 
features at u = Hl(r) ± nuJ o as shown in FigHJ^b). Besides, if the localized level is 
broadened so that there is a finite density of states around I^l(r), the electron can also 
bypass the SMT via the phonon-assisted resonant tunneling processes. According to 
our previous studies, [21] these phonon-assisted resonant tunneling processes can result 
in the phonon sidebands at uo = e =t n ^o, which also exhibit some sharp asymmetric 
lineshapes at uj = ^l{r) ± tiuq under extremely low temperature. These structures 
superpose over the Kondo-phonon-satellites, jointly responsible for the apparent sharp 
asymmetric peaks observed in Fig^b). 

The nonlinear differential conductance versus the bias voltage, i.e., G(Hws), is 
plotted in FigOl for very low temperature, where G reaches a peak value at V^ as = 
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and then decreases rapidly as Vuas increases. This zero-bias anomaly is usually regarded 
as the experimental evidence for the Kondo effect in the quantum dot systems. When 
eVbias = nhuo, there are some Kondo-phonon-satellite peaks show up, as evidenced 
by the very recent experiments. [15-17] According to the general understanding, in 
the non-EPI case the bias voltage plays a role similar to temperature, which means the 
differential conductance decreases logarithmically as Vb ias increases, [51] one may wonder 
here why the apparent sharp asymmetric features of the Phonon-Kondo-peaks in Fig|3] 
have not been smoothed out by bias voltage. It is because the Fermi distributions 
in the leads have sharp Fermi surfaces at very low temperature, the sharp features of 
the spectral function, as shown in Fig^b), can still be preserved even for finite bias 
voltage. As illustrated in the inset of FigEl once an asymmetric peak in the spectral 
function matches one of the Fermi levels, a sharp asymmetric differential conductance 
peak will appear. Therefore, the apparent sharp asymmetric features in GiVuas) are 
just the manifestations of the asymmetric lineshapes of the spectral function, which 
contain both contributions from the phonon-assisted Kondo processes and the phonon- 
assisted resonant tunneling processes as well. These sharp asymmetric features can only 
be smoothed out with the increasing temperature. 

It is worth pointing out the resemblances between our spin degenerate results and 
those by the real-time diagrammatic formulation [24, 25] or the generalized Schrieffer- 
Wolff transformation [30], where the discontiguous and asymmetric lineshape features 
of the Kondo-phonon-satellites in the spectral function or the differential conductance 
can also be observed, but have not yet been clearly pointed out and properly explained 
previously. 

4-2. finite Zeeman splitting case 

Now let us take a look at the spin non-degenerate case, where a Zeeman splitting of A 
for the localized level may be induced by a local magnetic field, i.e. = e + A/2 and 
e i = £o - A/2. 

Without the EPI, our calculation for the spin-splitting case shows that, compared 
to the spin- degeneracy case, the Kondo peak in the spin-resolved spectral function A^(oj) 
(or Ai(u)) generally decreases in magnitude and shifts away from fi by A (—A), which 
agrees with previous researches. [24-29,42,43] When the EPI turns on, as shown in 
FigEfa), the phonon satellite structure may develop in a distinct way, that is, for 
moderate Zeeman splitting only the satellites above (below) the main peak appear in 
A r (u) {Ai(u)). 

This peculiar feature in the spin-resolved spectral function stems from the interplay 
between the EPI and Kondo effect in the presence of the finite spin-splitting. It also 
manifests itself in the differential conductance. As shown in Fig0Jb), there is no satellite 
appearing in the region fi — A < eV^ as < /i + A, which is generally true as long as the 
system is at very low temperature. According to our simple picture above (c/FigJ2J), it 
is easy to understand. 
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Figure 4. (color online) The calculated total (blue dash dot line) and spin-resolved 
(black solid line for spin up, and red dash line for spin down) spectral functions in 
the equilibrium (a), and the nonlinear differential conductance (b) in the presence of 
finite Zeeman splitting. Here e = -2.5T, g = 0.6, k B T = 2 x l(T 4 r, U = 100r, 
?uj = O.ir, and A = 0.2I\ 
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Figure 5. (color online) The spin-resolved spectral function (black solid line for spin 
up, red dash line for spin down component) for different Zeeman splitting: (I) A = 0, 
(II) A = O.Oir, (III) A = 0.08I\ and (IV) A = 0.2T. Other parameters are taken as 
V blas = 0, e = -2.5r, g = 0.6, U = 100r, and haj = OAT. 



The spin-resolved spectral functions for four different Zeeman splitting values are 
plotted in FigEJ demonstrating that when the Zeeman splitting turns on, the overall 
Kondo peaks in the spin down spectral function shift towards the low energy, while 
those for spin up make a blue shift. The Kondo satellites above (below) the main peak 
in the spin down (up) spectral function decrease significantly as the Zeeman splitting is 
increased, then disappear when the A exceeds the line-width of the Kondo peak which 
is proportional to the Kondo temperature ksTx- However, for larger splitting, the main 
Kondo peak itself will be much suppressed, so that this spin-dependence of the Kondo 
phonon satellites will not be resolved at all. 

To understand how this pattern develops, it is helpful to recall that the satellite 
on different sides of the main peak is associated with different spin exchange process 
as depicted in Fig|2j This peculiar pattern for finite Zeeman splitting indicates that 




Figure 6. (color online) Maps for the spin up (a) and spin down (b) spectral functions 
as functions of uj and eVuas- The sections are plotted in (c) for four typical bias 
voltages: (I) eV bias = 0, (II) eV bias = 0.03I\ (III) eV Uas = 0.1T, (IV) eV bias = 0A5T. 
Here A = 0.2r, e = -2.5I\ g = 0.6, U = 100r, and huj a = 0.1T. 



the dominant contribution to the main peak in A^(u) comes from the coupling of the 
local state with the excited conduction electron states, while for A^(u) the dominant 
contribution comes from the coupling with the excited hole states. If A = 0, as 
illustrated in FigEl both spin exchange processes contribute to the Kondo main peak. 
When A starts to increase, both spin exchange processes decay in magnitude, but at 
different rates. When Zeeman splitting increases further, it would be possible that while 
one spin exchange process remains finite, the other is already totally suppressed. This 
results in the spin separation of the Kondo phonon satellites in the spectral function. 
However, when A k B T K , so that both processes are unlikely to happen and finally 
the Kondo peaks are quenched. 

The spin down (up) spectral function for fixed Zeeman splitting is shown in FigEfa) 
(FigEJb)) as function of u and Vu a s, i n which the highlighting cross lines correspond to 
the Kondo satellites. For explicitly, the sections at the four typical bias voltages marked 
by the dash lines in FigEfa) and FigJHfb) are also plotted in FigEJc), respectively. One 
can clearly see that, when the bias voltage turns on, each Kondo peak is separated into 
two sub-peaks related to two Fermi levels, respectively. These two sub-sets of Kondo 
satellites separate from each other further when Vu as increases, resulting in the peculiar 
pattern as shown in FigE^c). 



13 



5. Conclusion 

In summary, taking advantage of the improved treatment of the electron-phonon 
interaction and the extended equation of motion approach, we have systematically 
investigated the Kondo effect in the SMT in the presence of the electron-phonon 
interaction, finite bias voltage, and Zeeman splitting. The phonon-Kondo satellites 
observed in the recent SMT transport experiments [15-17] have been confirmed and 
explained by our results. Moreover, peculiar patterns of phonon-Kondo satellites in the 
spectral function as well as the differential conductance have also been predicted and 
explained by a clear physics picture. Although it is the interplay between the EPI and 
Kondo processes that discussed here, the physics picture and the theoretical calculation 
presented here may also be applicable to some other cases. For instance, the interplay 
of electronic internal excitations and the Kondo processes can also manifest themselves 
in some parallel Kondo satellites. [8] 
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